Local Time-Reversal Symmetry Breaking in d x 2_ y 2 Superconductors 
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We show that an isolated impurity in a spin singlet d x 2_ y 2 superconductor generates a d xy 
order parameter with locally broken time-reversal symmetry. The origin of this effect is a coupling 
between the d x 2_ y 2 and the d xy order parameter induced by spin-orbit scattering off the impurity. 
The signature of locally broken time-reversal symmetry is an induced orbital charge current near 
the impurity, which generates a localized magnetic field in the vicinity of the impurity. We present 
a microscopic theory for the impurity induced d xy component, discuss its spatial structure as well 
as the pattern of induced current and local magnetic field near the localized impurity spin. 
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There is now strong evidence to support the identifi- 
cation of the superconducting state of many of the high 
T c cuprates with a spin-singlet pairing amplitude hav- 
ing "d-wave" orbital symmetry, or more precisely d x 2_ y 2 
symmetry!!! This phase preserves time-reversal (T) sym- 
metry, but changes sign under reflection along the [110] 
and [110] mirror planes, as well as 7r/2-rotations in a 
tetragonal crystal. As a consequence d x i_ y i pairing cor- 
relations are particularly sensitive to scattering of quasi- 
particles on the Fermi surface. In this article we show 
that an isolated impurity in a spin singlet d x 2_ y 2 su- 
perconductor generates a complex d xy order parameter 
(OP) with locally broken T symmetry; the signature of 
this effect is an induced orbital charge current near the 
impurity and a localized magnetic field in the vicinity of 
the impurity. 

Atomic scale impurities, or defects, scatter conduction 
electrons which leads to local suppression of the super- 
conducting OP (pair-breaking) near the impurity. The 
mechanism responsible for pair-breaking is the forma- 
tion of quasiparticle states near the Fermi level which are 
bound to the impurity by Andreev scattering.0 The corre- 
sponding reduction in the spectral weight of the pair con- 
densate is responsible for pair-breaking. The existence 
of quasiparticle states near the Fermi level can also lead 
to local Fermi-surface instabilities and mixing of order 
parameters with different symmetry.! 3 ! Low temperature 
phase transitions associated with a secondary OP may 
provide new information on the mechanism(s) for pair- 
ing in unconventional superconductors ,□ while impurity- 
induced mixing of the d x 2_ y 2 OP with an OP of different 
symmetry can provide direct information on the atomic 
structure of the impurity.B 

Recent transport experiments report evidence for low 
temperature phases associated with a secondary OP at 
surfaces (T s ~ 8K).Q This was interpreted in terms 
of a surface phase transition to a cLri-j,2 + is state 
with spontaneously broken T-symmetry.a Bulk transport 
measurements on Bi2Sr2CaCu20s+5 (Bi-2212) show a 



LA-UR: 99-1349 



pronounced drop in the thermal conductivity at_T* ~ 
150 mK < T c w 80 K in Ni doped Bi-2212.Q This 
anomaly was interpreted as the signature of a second 
superconducting phase with a fully gapped quasiparti- 
cle spectrum and a mixed symmetry OP of the form 



id r 



This phase was proposed to arise from 



a coupling of the orbital momentum of the conduc- 
tion electron with the spin of the magnetic impurity, 
7~L S o = J dr^l t (r) v(r) L orbit • S lmp ^Q,(r).H Measurements 
of the spin-orbit coupling energy for Ni 2+ ions indicate 
that it is a few n percent of the nonmagnetic and ex- 
change channels.Ea In this model the id xy OP is induced 
at a temperature above the second phase transition, 
T* < T <T C ; the low temperature transition is argued to 
be ordering of the impurity-induced "patches" of the lo- 
cal d x 2_ y 2 Hd xy order Above T* patches with randomly 
fluctuating internal phase destroy the long range order, 
(d X y) = 0, but preserve {\d xy \ 2 ) ^ 0. Thus, the local 
structure associated with d x 2_ y 2 ±id xy symmetry near a 
magnetic impurity should be observable at temperatures 
well above T*. The electronic and magnetic structure 
near an impurity located near the surface of a supercon- 
ductor can now be studied with atomic resolution at low 
temperature by scanning tunneling microscopy (STM), 
opening a new window for local probes. BIMj 

In this article we investigate theoretically the local 
structure of the OP and the current distribution in the 
neighborhood of an atomic impurity within the d x 2_ y 2 
model for the high T c cuprates. We present new ana- 
lytical and numerical results for the local d x 2_ y 2 ± id xy 
OP that develops near a magnetic impurity, which has 
attracted new theoretical interests. OIL 2 ! Our approach 
follows closely the theory developed in the late 70's 
for ions in superfluid 3 He,E 3 J!3 and later adapted to 
study the properties of impurities in heavy fermion 
superconductors.Oc3 The theory of impurity scattering 
in superconductors can be formulated to quasiclassical 
accuracy as an expansion in ct/£oj where a is the (lin- 
ear in 2D) cross-section of the impurity for scattering 
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of normal-state quasiparticles at the Fermi surface, and 
Co = Vf/irAo is the superconducting coherence length. 
This ratio is typically very small in low T c superconduc- 
tors and in superfluid 3 He, but may be as big as 0.2 for 
strong scatterers in high T c superconductors. 

We start from Eilenberger's transport equation, 
for the matrix propagator in particle- hole/spin space, 
g(Pf, r; The diagonal element of the propagator de- 
termines the local density of states and local equilibrium 
current distribution, and the off-diagonal elements are 
the components of the local pair amplitude. Quasiparti- 
cle scattering off an isolated impurity is included throug 
a source term on the r.h.s of the Eilenberger equation 
The transport equation can be linearized to leading order 
in er/£o for distances r > <r from the impurity. In this 
limit the Fourier transform of the linearized transport 
equation reduces to,E3 



ie n T 3 -A b -a imp ,6g +q-v f 5g= i+SA + Sa 



, 9b 



(1) 



The bulk propagator, <7&(p/;e n ) = — 7r[ie n f 3 - A b (j>f)]/E, 
order parameter, Ab(p/) = Ab(p_f)fi ia 2 , impurity scat- 
tering self-energy, a lmp , and in-plane Fermi velocity, v/, 
are inputs to the linear response equations. The energy 
denominator is given by E = (|Af,(p/)| 2 + e 2 ) 1 / 2 , where 
+ |Trf3CTi mp (e n ) is the renormalized Matsubara 
frequency. The t matrix for the isolated impurity, as well 
as the induced OP, <5A = [SA^i + SA 2 f 2 ]ia 2 , an d self- 
energy, 5<7i mp , enter the r.h.s. of Eq. (|l|) as source terms. 
Here fj and Oi are Pauli matrices in particle-hole and 
spin space, respectively. The t matrix for the isolated 
impurity is given by 

i(p f ,p' f ;e n ) =v(p f ,p' f )+Nf J d 2 p'}v(p f , p" f ) 

X9b(Pf;e n )t(p'},p' f ;e n ), (2) 

where Nf is the 2D density of states at the Fermi energy 
per spin, and u(p/,pj) is the impurity potential, which 
is evaluated in the forward scattering limit in Eq. (jl]); 
p'j = pf. We separate v into channels for nonmagnetic 
(u), spin-spin exchange (m = JSj mp ), and spin-orbit 
scattering (u so ) between the orbital momentum of the 
quasiparticle and the impurity spin, S; mp , as well as the 
self-coupling (w so ) to the quasiparticle spin, 

w(p/,P/) = «(p/,P/) + m(p/,p}) • S + 

*[ u so(p/,P/) • S imp + w so (p / ,p' / ) • S]f 3 . (3) 

The induced OP is self-consistently determined from the 
gap equation, 

6A(p f ,q)=N f T^2 JcPp' f V(p f ,p' f )f(p' f ,q;e n ), (4) 



V 2 rj 2 ((j))ri 2 ((f)'), where the eigenfunctions are ?7i(</>) = 
cos2</> and rj 2 ((f) — sm2cf> for the two channels, respec- 
tively. The dominant, attractive interaction is V± = 
V x 2_y2, and the subdominant interaction, V 2 = V xy , may 
be either attractive or repulsive. We neglect the s-wave 
pairing channel in order to simplify the analysis, and we 
restrict our discussion to the regime in which the sub- 
dominant interaction, V xy , is repulsive or too weak to 
nucleate a bulk d xy OP. 

Numerical calculations of the OP and current dis- 
tribution were carried out for an isolated impurity at 
r = 0, with the t matrix source term of the form 
t(pf, p/; e n )S(r) in real space. We modeled the posi- 
tion of the impurity to quasiclassical accuracy by replac- 
ing the delta function by a smooth function, 5 ro (r) 

- exp(— r 2 /tq), of atomic width 7'o = 0.1£o. This model 
guarantees a smooth cutoff in q-space and faster con- 
vergence of the Fourier integrals. For the computation 
reported here we also chose a subdominant pairing inter- 
action corresponding to a bare subdominant transition 
temperature of T c2 /T c \ — 0.1, which is well below the 
threshold for bulk nucleation of a d xy order parameter. 

The physical solution to the linearized transport equa- 
tion (0) isjia 



^(p/jq; £r. 



E g b + ttQ 
2tt(E 2 + Q 2 ) 



[t + SA + Sa imp , g b ] , (5) 



where Q — iq • v/. The induced charge current is also 
determined by the t matrix and induced OP, 



Sj(<D=N f Tj2 d 2 p f v f 



2z7reQA 1 
E(E 2 + Q 2 ) 



(t 2 +6A 2 ), (6) 



where t 2 is the f 2 component of i. 

We evaluate the t matrix in second-order Born approx- 
imation, which is adequate for weak scattering. More 
importantly, the Born approximation generates the rel- 
evant coupling between the d x 2_ y 2 and d xy order pa- 
rameters. We also assume that the impurity poten- 
tial is short-ranged, so we retain only the scattering 
amplitudes for the s-wave and p-wave channels, i.e., 
u(Pf,p' f ) ~ u + mpf ■ p' f , and m(p / ,p' / ) w (J + 
J\Pf ■ p'y)Sj mp . For the spin-orbit terms we approx- 
imate, u so (p/,p^) w (A + Aip/ ■ p' f )pf x p' f , and 
w so(p/,P/) ~ (wo+ w iP/'P/)P/ X P/- The £ matrix then 
has the general form i = \tirjiTi-\-t 2 rj 2 T 2 \ i ia 2 + [t3+t-S] f 3 . 
The important interference term is given by 



dcj> A ir] 2 (4>) 
27r v /|A m (0)| 2 +e~ 2 ' 



(7) 



where V(p/,Py) is the pairing interaction and / = 
[<5/]Ti + Sf 2 r 2 ]ia 2 is the induced off-diagonal pair am- 
plitude. We resolve the pairing interaction into the dom- 
inant attractive d x 2_ y 2 channel, and a secondary pairing 
channel with d xy symmetry, V(pf,p'f) = Vir/i (^)ryi {4>') + 



with the spin-orbit parameter A 2 = itoAi +""iAo + Jo^i + 
J1W0, and the impurity spin S z . The t 2 term generates 
a correction to the off-diagonal propagator with d x 2_ y 2 
(Big) symmetry and induces the id xy (B 2s ) OP near the 
impurity. 
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FIG. 1. Contour plots of the induced id xy OP (top left), 
the current density (top right), and the magnetic field (bot- 
tom left), at T = 0.1T C for T c2 /T cl = 0.1 with no bulk scat- 
terers. Solid (dashed) lines are negative (positive) contour 
lines. Bottom right: Field plot of the current density. 
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FIG. 2. Induced id xy OP at q = 
by the Ai(T) OP vs. T c 
ing rates a = 1/(2ttT c it] 
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(a) SA22 normalized 
= T c (a) for different bulk scatter- 
(b) SA22 and the bulk A 2 OP for 



attractive pairing potentials (T C 2/T c i) at T = for a = 0.0 
(solid) and a = 0.1 (dashed). See text for the definition of S. 

The corrections to the bulk OP, <5A,, with i = 1,2, do 
not belong to a single representation, i.e., <5Ai(c/>,q) 7^ 
r]i(4>). In particular, <5A 2 (</>,q) has mixed and d x 2_ y 2 
symmetry, <5A 2 (0,q) = <JA 2 i(q)7?i(0) + (JA 22 (q)?72 (</>)■ 
These amplitudes also determine the current distribution 
(H), and satisfy the coupled equations 



(-)^A 2l ( 9l , -g 2 ) = (-ySA 2l (q 2l qi ). The induced <SA 22 
OP component with B 2g symmetry is finite for q = 0, 
while the induced 6A 2 i component with Bi g symmetry 
vanishes for q = and along the diagonals and princi- 
ple axes. The Fourier transformation of <5A 2i (q) also 
vanishes, i.e., <5A 2 i(r) = 0; thus only the induced compo- 
nent with B 2g symmetry survives. The contour plot of 
5A 2 (r) in Fig. 1 shows a four- fold pattern characteristic 
of the d xy amplitude with maxima located at approx- 
imately 0.3£o along the nodal directions of the d x 2_ y 2 
OP. 

Bulk impurity scattering is pair-breaking for any un- 
conventional OP including the induced d xy amplitude. 
Figure [2] shows both the temperature dependence of the 
induced d xy OP and the pair-breaking suppression of 
<5A 22 (q = 0) by bulk impurity scattering. Note that 
the d xy OP develops below T c and that it is suppressed 
by bulk scattering on the same scale as the bulk d x 2_ y 2 
OP. Fig. |(b) shows the increase in the induced d xv 
amplitude with increasing (attractive) pairing interac- 
tion in the B 2g channel (T c2 /T c i); the divergence at 
TcijT c \ > 0.37 corresponds to a bulk instability for 
d x 2_ y 2 ± id xy pairing. For a repulsive pairing interaction 
V 2 < we find a cutoff dependent result for the induced 
OP, -0.44££g/log(w c /Ai) < £A 22 (g = 0) < 0, at T = 0. 
We parameterized the interference term in the scattering 
amplitude by the coupling energy £ — -kS z Nj\ 2 /£g, 

The existence of a Hd xy OP implies that the equilib- 
rium superconducting state breaks T symmetry locally 
near the impurity, in addition to broken [110] and [100] 
reflection symmetries. The signature of the d x 2_ y 2 ±id xy 
state is the equilibrium charge .cur-cent and magnetic field 
distribution near the impurity.E3iiju3 From the t matrix 
in Eq. (|J) and the induced d xy OP we obtain 
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FIG. 3. Spatial dependence of the normalized OP, 
<5Aj ; (r), current density, <5j' n '(r) (x and y components), 
and magnetic field, Bj 1 (r), along the (010) and (110) di- 
rections at T = 0.1T C for a modest d xy pairing interaction of 
T c2 /T c i — 0.1 (well below the bulk transition to d x 2_ y 2+id xy ). 



[1/ Vl - £n (q)] <5A 21 (q) - JC 12 S A 22 (q) = A 1 (q) 
[I/V2 - £22 (q)] <5A 22 (q) - /C 12 SA 21 (q) = A 2 {q) 



(8) 
(9) 

and 
The 



where /C y (q) - ttT^ J + Q : 

A(q) = irT E £ „ h(e n ) J p^^E / [E 2 + Q 2 
solutions to these equations have, in addition to in 
version symmetry, the mirror reflections SA 2 i(qi,q 2 ) = 



6j(ci)=iireN f Tj2 J gv/(» v/(0)-q- 



Ai??i(0) 



l E(E 2 + Q 2 ) 

x Ut 2 (e n ) + ^A 22 (q))?7 2 (0) + 6A 21 (q) m (0)) . (10) 



It is straightforward to show that the current den- 
sity obeys the symmetry relations: <5j(q) = — Sj(— q), 
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and Sj l (q 1 ,q 2 ) = {-)' l+1 5j. l {-q 1 ,q 2 ), Sj 1 (q 1 ,q 2 ) = 
~5j 2 {q 2l qi). From these relations one might expect a 
simple circulation pattern for the induced charge cur- 
rent; however, as Figs. 1 and || show, the current den- 
sity exhibits the superposition of a very small circulat- 
ing current loop and four counter circulating currents 
around the nodal directions, which are anchored to the 
local maxima of the induced SA 2 (r) OP. This pattern is 
qualitatively similar to the current-distribution predicted 
for a vortex with d+is symmetry,L3 however, there is no 
circulation at large distances from the impurity. The 
complex flow pattern that is observed near the impurity 
is also observed for mesoscopic superconductors with sur- 
faces that are normal to the (110) direction,E3 and reflects 
the strong nonlocality of the current response shown in 
Eq. ©. 

The spatial pattern of current generates a four-fold 
magnetic field distribution which we calculated from the 
current distribution using the Biot-Savart law, B z {r) = 
\ J d 2 r'\r' - r|- 3 (r' - r) x <5j(r'). The field distribution 
in Fig. 1 shows 8 sectors of flux threading in and out of 
the plane. As a result the net magnetic flux through the 
superconducting plane is zero. This fact was checked nu- 
merically. The magnetic flux was calculated for squares 
of area L 2 and shown to vanish in the limit L — > 00. 
This is a general result, at least at the quasiclassical level, 
provided that scattering by the impurity does not gen- 
erate-, a coupling of the current to a soft mode of the 
OP £-3 However, particle-hole asymmetry corrections to 
the quasiclassical theory may lead to a net moment from 
the impurity-induced orbital currents. 

The magnitudes of the induced id xy OP and magnetic 
field near an impurity depend on parameters characteriz- 
ing the interaction between quasiparticles and the mag- 
netic impurity. Not much is known about these interac- 
tions. Thus, measurements of the induced OP or mag- 
netic field near an impurity can provide direct informa- 
tion about the coupling of the quasiparticle orbital mo- 
mentum to the impurity moment. We can express the 
impurity induced OP, current density and magnetic field 
in terms of a few key material parameters of the impu- 
rity. We scale the induced OP in units of the coupling 
energy, 5A { 2 n \r) = SA 2 (r)/£. From Eq. @ we obtain 
the scale of the current and field: <5j(")(r) = <5j(r)/ (cB), 



BT'(v) = B z (r)/B, with B = ZN f v f £ = jf^^j. Fig- 
ure H shows the spatial variations of the scaled OP, cur- 
rent density and field profile along the (110) and (010) 
directions at low temperature. 

We estimate the magnitude of the induced OP and 
magnetic field for Ni impurities in Bi-2212 as follows: 
The coupling parameter for Ni 2+ ions is estimated from 
the spin-orbit coupling energy for free Ni 2+ ions,B A 2 ~ 
(30meV • a 2 ,) 2 . The in-plane penetration depth of A ss 
200 nm, the interlayer spacing, d 1.5 nm, the in-plane 
lattice constant, ao ~ 5.4 A, and the Fermi velocity, 
Vf ~ lOOkm/s, provide a determination of the density 
of states per Cu-0 bilayer, Nf = c 2 d/(Aire 2 v 2 X 2 ). This 
gives an estimated energy scale for the induced d xy gap 
of £ ~ 0.1 mcV, and a magnetic field scale of order 



y 2 ± id xy phase order transition at lower tempera- 
K, as observed in 0.6% and 1.5% Ni 



B ~ 1 fiT. An induced d xy gap of order 0.1 meV w 
IK, is the right order of magnitude to account for a 

d 2 

ture, T* ~ 15 
doped Bi-2212. 

In conclusion, we have shown that spin-orbit scatter- 
ing induces a d x i_ y i +id xy state, which locally breaks T 
symmetry in the vicinity of a magnetic impurity. The in- 
duced OP develops below T c and survives bulk impurity 
scattering so long as the bulk d x 2_ y 2 OP does. The signa- 
ture of the spontaneously broken T symmetry manifests 
itself as a complex pattern of circulating charge currents 
near the local maxima of the d xy OP located along the 
(110) and (110) directions. We estimated the magnitude 
of the induced d xy gap to be ~ 1 K, which should be ob- 
servable in low temperature STM measurements of the 
tunneling density of states. 
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